Poincare recurrences from the perspective of transient chaos 
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We obtain a description of the Poincare recurrences of chaotic systems in terms of the ergodic 
theory of transient chaos. It is based on the equivalence between the recurrence time distribution 
and an escape time distribution obtained by leaking the system and taking a special initial ensemble. 
This ensemble is atypical in terms of the natural measure of the leaked system, the conditionally 
invariant measure. Accordingly, for general initial ensembles, the average recurrence and escape 
times are different. However, we show that the decay rate of these distributions is always the same. 
Our results remain valid for Hamiltonian systems with mixed phase space and validate a split of 
the chaotic saddle in hyperbolic and non-hyperbolic components. 
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The idea of recurrence to a certain region of the phase 
space was introduced by Poincare in his famous studies 
of the three body problem as a stabiUty criterion [l[ . The 
wider implications of his results were soon recognized and 
Poincare recurrences have played an important role ever 
since the first debates on the foundation of nonequilib- 
rium processes 0, 0, 01 • More recently, recurrences have 
become a standard tool to investi gate low-dimensional 
closed chaotic systems d 0, 0, i, i,M El, [ll [ll ■ 

Open systems, on the other hand, are usually inves- 
tigated as scattering problems: the long-lasting chaotic 
transients are related to an invariant fractal saddle 
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16[ and only short time escapes depend on the initial 
ensemble. A sharp distinction between open and closed 
systems exists since the time of Poincare [l| . 

Closed systems can be converted into open ones by 
defining a finite region of the phase space as a leak. 
Leaking dynamical systems mimics the effect of exper- 
imental observations 17 , [l^ US] and has also been ap- 
plied as a tool to investigate the dynamics of closed sys- 
tems |2C 
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In this Letter, we find the complete correspondence 
between the recurrence and the leak problems. We then 
apply the transient chaos theory of open systems to de- 
scribe Poincare recurrences of closed systems. More pre- 
cisely, considering the recurrence region to be the leak, 
we show that for chaotic (both dissipative and Hamilto- 
nian) systems (i) the exponential relaxation rate of the 
recurrence problem and the escape rate from the leaked 
system are always identical, (ii) the relaxation rate can 
be expressed by the conditionally invariant measure Q 
of the leaked region, (iii) under a properly chosen ini- 
tial ensemble of the escape problem, the entire recur- 
rence and escape time distributions coincide, and (iv) 
the points above remain valid even in Hamiltonian sys- 
tems with mixed phase space, where for recurrence/leak 
regions far from Kolmogorov-Arnold-Moser (KAM) is- 
lands 15, [3] the escape rate and an associated hyper- 



bolic saddle are shown to be well defined for times shorter 
than a crossover time. 

To be specific, consider a chaotic map M defined on a 
bounded phase space P with an invariant chaotic set: a 
chaotic attractor or a chaotic sea. The natural ergodic 
measure on this invariant set is denoted by its den- 
sity by Pfj^. The recurrence region is a subset / C P. 
We do not restrict ourselves to small ^(I) [13, 11, [13] ■ 
The Poincare recurrence theorem 0, 0, 0| ensures that 
for almost any initial condition in /, there is a first re- 
currence at some discrete time ni, a second one at time 
TT-i -\- 712, etc. In the long time limit, the set of recur- 
rence times Hi defines the recurrence time distribution 
Pr{n). For chaotic systems 0,13, El, this distribution is 
exponential for n larger than some small n* 
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Prefactor g^ and the relaxation rate 7^ depend on the 
choice of /. While the average recurrence time {n)r is 
given by Kac's lemma [3]: 



{n)r 
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it has been unknown whether the relaxation rate 7^ can 
be expressed by means of any measure of /. The equality 
7r — 1/ {ur) = pXI) is proved only in the unrealistic limit 
fi{I) ^ (see [sl, and references therein). 

Open (transiently) chaotic maps are characterized by 
the escape time distribution Pe{n), which gives the frac- 
tion of trajectories of a certain initial ensemble that 
leaves the system at time n. This distribution is also 
exponential [l^, [l6|. We are interested in the class of 
open systems obtained by leaking a closed system at a re- 
gion / C P, described by map M''{x) = M{x) for x G T\I 
and M^x) ='exit' for a; e /. Notice that the escape hap- 
pens one step after entering / and initial conditions can 
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be in /. In particular, for such a leaked system we have 

PeW«5ee-^=", (3) 

for n > n*, where 7e is the escape rate. The average 
escape time (n)e, also called the lifetime of chaos, is usu- 
ally estimated by the reciprocal of 7e 15|, m. While 
the full Pe{n), and the average {n)e, depend also on the 
choice of the initial density po used in generating Pe(n), 
the escape rate does not (provided it exists). All these 
quantities depend on /. The theory of transient chaos 
isl . 3] explains these results based on the existence of 
an invariant chaotic saddle in the complement set of the 
leak T\I The saddle is the union of all trajectories 
that never enter the leak, neither forward nor backward 
in time. It is a fractal set of measure zero for strongly 
chaotic systems. Trajectories spending a long time out- 
side / must come close to the saddle along its stable man- 
ifold, stay in the vicinity of the saddle, and leave it along 
its unstable manifold. Therefore, irrespective of po, the 
long-term emptying process is governed by the invariant 
saddle, which is uniquely defined by /. 

We now apply this picture to Poincare recurrences. 
Long recurrences should correspond to trajectories that 
fall near to the chaotic saddle's stable manifold right af- 
ter exiting the recurrence region /. Therefore, we expect 
the relaxation rate of the recurrence distributions 7^ and 
the decay rate of the escape distribution 7e to coincide 



7r = 7e 



7, 
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for a given recurrence/leak region /. This equality is 
illustrated for the Henon map in Fig. [U where it is ap- 
parent that the slopes of liipr{n) vs. n and Inpe(ri) vs. n 
coincide. Equation (j4]), which will be rigorously justified 
below, implies that to any Poincare recurrence problem 
there exists a chaotic saddle (that of the corresponding 
leaked system) that governs the long-term recurrences 
and determines 7^. This saddle is shown in Fig. [^b- Its 
double fractal character makes evident the improvement 
of Eq. (Ill) in comparison to the previous results being 
limited to p{I) ^ 

Despite having the same exponential decay (|4]), Pr{n) 
and Pe{n) are different. This holds for practically any 
initial density pQ. We show now that there is, neverthe- 
less, a special pa = pr 0, of the leaked system for 
which pe(n) —pr{n). Consider inserting trajectories into 
the system through the leak. More precisely, consider 
the density pr{x) obtained as the first iterate through 
map M of the points x £ I distributed according to the 
natural density of the closed system p^: 

Prix) = p^{M-\x) n I)/p{M{I)) for X e M{I), (5) 

where M^^{x) fl / denotes the points that come from / 
and p{M{I)) ensures normalization. This distribution 
is shown in Fig. [5^. In view of the ergodic theorem. 




FIG. 1: (Color online) (a) Recurrence Pr{n) (B/above) and 
escape Pe{n) (A/below) time distributions in the Henon map 
shown in Fig. (2] For pr{n) a single initial condition (xq, yo) = 
{xc,yc) was iterated over 10^^ steps. For Pe{n), po ~ Pi_i 
was built by 7 • lO'' trajectories. Inset: short time behavior, 
where Pe{n) with po = Pr, given by ([5]), is marked by O. 
(b) pe,rexp(7n), with 7 = 0.055 [p4I) = 0.0535 < p(l) = 
0.0377]. Numerical limitations appear for large n. 




FIG. 2: (Color online) Phase space of the Henon map 
{x„+i,yn+i) = (1 — lAx^ + y„,0.3x„) with 7 as a circle of 
radius 0.05 centered at (sc, yc) = (1.0804,0.0916). (a) Henon 
attractor: support of p^j (thin line) and of pr (•) at M{I), 
see ([5}. Inset: pr{x,y) projected onto x. (b) Chaotic saddle 
(x) and the support of pc (thin line). 



the distribution of all positions of a single infinitely long 
recurrent trajectory, one iteration after returning to /, is 
precisely given by pr- Since the escape times for pr are 
exactly the recurrence times, it follows that 

Pr{n) ^Pe{n) with pa = Pr, 

for all times n (see the inset of Fig. [T^). Of course 7^ — 7e 
and, because 7e is independent of po, this implies 

Having shown how recurrences can be viewed as es- 
capes from a leaked system, we now perform a formal de- 
scription in terms of the ergodic theory of transient chaos. 
A central concept in this theory is the so-called condition- 
ally invariant measure, c- measure in brief [2^ . 25 . [2^ . 
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Smooth initial densities iterated via the action of map 
do not converge to any finite measure due to the 
permanent escape. If, however, one applies a compensa- 
tion in the form of multiplying the density in each step 
by exp (7) > 1 it converges to a finite measure, the c- 
measure /ic- Due to the permanent contraction along 
the stable direction, the density pc of the c-measure is 
nonzero along the unstable manifold of the saddle [2^. 
Qualitatively speaking, this measure characterizes the es- 
caping process from the saddle and is thus completely 
different from the natural measure of the closed system. 
With the normalization /ic(r) = the compensation 
factor exp(7) is the c-measure of the region not escaping 
within one time step [23I]. In the leaked system this is 
^c(r \ /) = A^c(r) ~ ^J'c{I) — ^ — tJ-cil) and therefore 

7 = - In (1 - p,{I)) [« for small (6) 

Note that the c-measure itself depends on the form, size, 
and location of /. In view of (|4]), it is remarkable that 
the relaxation rate 7^ of the Poincare recurrences of the 
closed system is given by the c-measure of the recurrence 
region, when viewed as a leaked system. Furthermore, 
Po = Pc is the only initial density for which the exponen- 
tial decay starts from the very first iterate (n* = 1). As 
a consequence, 

(n)e = l/pc{I) with po = Pc- 

The reason for the strong difference between this rela- 
tion and Kac's lemma ^ (valid for recurrences and for 
Po — Pr) is that pr in ([5]) is atypical from the point of 
view of the c-measure. In fact, pc is concentrated along 
the saddle's unstable manifold while all points on pr 
have come from the leak. Typical initial densities po 
quickly converge to pc, explaining why typical (tVL-s, but 
not {n)r, are well estimated by I/7 ~ 1/ pc{I) [27|. 

We illustrate now through two examples the new per- 
spectives opened by our unified approach, (i) Recently, 
it has been surprisingly observed that different leaks / 

Location 
In 
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with equal p{I) lead to different 7e's 
dependence of 7^ has also been reported 
view of (|4]), these works refer to the same phenomenon. 
Furthermore, Eq. ^ explains the observations by stat- 
ing that it is the c-measure and not the natural measure 
that determines 7. (ii) A basic problem in the mathemat- 
ical approach to recurrences is the convergence of Pr (n) 
to a Poissonian [7 = p{I)] for p{I) — > [9!]. A funda- 
mental property of the c-measure is its converges to the 
natural measure [i^ [2^ . Equation ^ implies that both 
convergences occur precisely in the same way. 

We turn now to Hamiltonian systems with mixed phase 
space, where deviations from the exponential decay ap- 
pear for long times. The trajectories that approach the 
KAM islands (or other non-hyperbolic structures) stick 
to them for a time that is roughly algebraically dis- 
tributed d, d, [l3|. This phenomenon is usually inves- 
tigated using recurrences 1, [l3| but appears also in 



escape problems [16|, l28|. The asymptotic power-law like 
decay of Pr^ein) is the same and independent of the choice 
of initial densities (away from islands [11]), because it 
is related to the non-hyperbolic structures. For shorter 
times, an exponential decay of Pr^e{n) is observed 
what has deserved little attention from the recurrence 
perspective so far. 



^ escape 
recurrence 




FIG. 3: (Color online) (a) Recurrence Priji) (B/below) and 
escape pe{n) (A/above) time distributions in the standard 
map shown in Fig. |4]with e — 0.1. For pr{n) a single initial 
condition was iterated over 10^^ steps. For Pe{n), po = P/j 
at |a:;| > 0.25 was built by 2 • 10^ trajectories. Lower inset: 
short time behavior. Upper inset: scaling of nc7 with 7, 
obtained by changing e in 0.03 < £ < 0.4. (b) pe.r exp{'yn) , 
with 7 = 0.0108 [p^I) = 0.0107 < p{I) = 0.0113]. 

Accordingly, for / outside KAM islands and for n > n* 

e^''" for n* < n < Ua 



PrA'^) 



ae -|- b{'-fn) " for n > Ua, 



(7) 



where ae~'^"° ^ b (7/10) The recurrence and escape 
problems can be again related using pr in However, 
the analysis based on the chaotic saddle has to be re- 
considered in the presence of KAM islands. Numerical 
results for the standard map in Fig. [3] indicate that, for 
intermediate times, jr,e are well defined and equal, as 
in Q. These results remain valid [27| for a wide range 
of Fs and of smooth po's, away from islands (see [2§| for 
po's touching the island). 

To understand these findings we propose here to effec- 
tively split the saddle in a hyperbolic (outside islands) 
and a nonhyperbolic (nearly space filling [l3|) compo- 
nent. To justify this splitting, it is worth defining a 
crossover time > n„ as ae~^"'' = 6(7nc)~", starting 
from which the nonhyperbolic component dominates. In 
the upper inset of Fig. [3l we present evidence for 



~ 1/7 [~ ^/ p{I) for small p{I)]- 



(8) 



This scaling can be obtained from the definition of ric by 
assuming a weak dependence of b/a on 7. Equation ([5]) 
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FIG. 4: (Color online) (a) Phase space of the standard map 
{xn+i,yn+i) = (xn + j/n+i.^n - 0.52 sin(27ra;„) ) , periodic in 
\x,y\ < 0.5. / = [0.25 < a; < 0.25 + e, -0.5 <y < -0.5 + e], 
e = 0.2 (ffl). The KAM island (non-hyperbolic component of 
the saddle) is at the center. The support of pr is at M{I) 
(dark/red region), (b) The hyperbolic component of the sad- 
dle generated by the method of [igI ]. using n — 60 < ric = 160. 



impUes that the exponential decay is always dominant for 
small / and that hyperbolic saddles and c-measures can 
be effectively defined for intermediate times n* < n < Uc 
(see Fig. \^). An important difference between the 
manifolds of the hyperbolic and nonhyperbolic compo- 
nents is that they attract/repel exponentially and sub- 
exponentially, respectively. The picture that emerges is 
that po relaxes to the non-hyperbolic component only af- 
ter, and mainly via, the hyperbolic one. Therefore, b/a 
hardly depends on po and, consequently, neither does Uc- 
This explains scaling ([8]) and also why the absolute value 
of Uc is approximately the same for escape and recurrence 
(see the upper inset of Fig. [S]). 

In summary, we have presented a new interpretation 
of recurrences based on the theory of transient chaos. In 
particular, we have expressed the exponential relaxation 
of the recurrence time distribution in terms of properties 
of chaotic saddles. In addition, in Hamiltonian systems 
with mixed phase space we could properly order the ex- 
ponential and power-law decays, respectively, to hyper- 
bolic and non-hyperbolic components of the saddle. Our 
results hold in continuous time and in any dimension, 
being thus valid, e.g., in classical many particle systems. 

So far we have emphasized that Poincare recurrences 
can be viewed as an escape problem. However, in view of 
the correspondence between the two problerns presented 
here, also the results on recurrences [1, 0, [13, El, 12 1 
apply to leaked systems. Furthermore, the split of the 
chaotic saddle in Hamiltonian systems with mixed phase 
space provides a theoretical description with practical 
applications in problems like resetting in hydrodynam- 
ics [iOl and emission from optical cavities 
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